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ABSTRACT 

One of the strongest predictions of the ACDM cosmological model is the presence of dark satellites orbiting 
all types of galaxies. We focus here on the dynamical effects of such satellites on disky dwarf galaxies, and 
demonstrate that these encounters can be dramatic. Although mergers with M sa t > Aid are not very common, 
because of the lower baryonic content they occur much more frequently on the dwarf scale than for L*-galaxies. 
As an example, we present a numerical simulation of a 20% (virial) mass ratio merger between a dark satellite and 
a disky dwarf (akin to the Fornax dwarf galaxy in luminosity) that shows that the merger remnant has a spheroidal 
morphology. We conclude that perturbations by dark satellites provide a plausible path for the formation of 
dSph systems and also could trigger starbursts in gas rich dwarf galaxies. Therefore the transition from disky 
to the often amorphous, irregular, or spheroidal morphologies of dwarfs could be a natural consequence of the 
dynamical heating of hitherto unobservable dark satellites. 

Subject headings: galaxies: dwarf, interactions, evolution; (cosmology:) dark matter 



1. INTRODUCTION 

According to the ACDM scenario, stellar disks are immersed 
in dark matter halos and are surrounded by a full spectrum of 
satellite companions. Encounters with these satellites can in- 
ject significant amounts of energy into the system, with con- 
sequences that vary from negligible to fully catastrophic disk 
destruction depending on the relative mass of the perturber and 
the configuration of the event (relative distances and velocities). 
Disk heati ng by such substructures has been addressed in previ- 
ous work (iToth & Ostrikeril 19921: lOuinn et al.lll993t iFontet alj 
1200 It iBenson et al.fl 2004). but has generally focused on the ef- 
fect on bright Milky Way-like galaxies. 

Cold dark matter models predict the structure of halos to be 
self-similar; in such a way that, when properly scaled, a Milky 
Way-si zed halo looks comparable to one hosting a faint dwarf 
galaxy dMoore et al.ll999tlSpringel et al.l2008tlKlimentowski et alj 
2010; IWang et alJl2012l) . However, galaxy formation is not a 
self-similar process, as the properties of galaxies depend in a 
complex way on e.g. the mass of their host halos. For ex- 
ample, low mass (dwarf) galax i es are much more inefficient 
at for ming stars telanton et al.1 120011 ; iRobertson & Kravtsovl 
2008) an d have much higher mass-to-light r atios than larger 
galaxies dYang et al.ll2003l:lWalker et alj|2009l ). In addition, gas 
cooling is likely to be (nearly) completely inhibit ed in dark mat- 
ter h alos with masses below ~ 10 8 /i _1 M (Kauf mann et al.l 
2007), which implies that the satellites of dwarfs should be gen- 
erally completely dark in contrast to satellites in galaxy clusters 
or around L*-galaxies. 

In this Letter we show that these considerations imply that 
the dynamical perturbations of dark-matter satellites on dwarf 
galaxies are much more important than on -galaxies. Dark 
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satellites may provide a channel for the formation of dwarf 
spheroi dal galaxies without the need to recur to en vironmental 
effect s dMaverll2010l) or multiple body interactions dSales et al.l 
2007). Such interactions may also be responsible for the ob- 
served increase of disk "thi ckness" towards fainter galaxies 
dYoachim & Dalcantonll2006i) , as well as explain the existence 
of isolated dw arfs undergoing intense starbursts without an ap- 
parent trigger (Ber gvalll20ll as a result of a major merger with 
a dark companion (T. K. Starkenburg et al., in prep.). 

2. MODELS 

Our goal is to quantify the effects of substructures on disk- 
like galaxies over a broader region of parameter space (and 
specifically mass range) than done in previous work. To this 
end, we use the second resol ution level of the Aquarius Sim- 
ulations (Springel et al. 2008) and study the assembly history 
of main (as opposed to satellite) dark matter halos in the mass 
range 10 8 - 1O 12 /i _1 M . We follow their evolution from t = 2 
Gyr onwards (z < 3), since by this time all halos in our sample 
have accreted at least 10% of their final mass, and the concept 
of "main/host" is well-defined. At the present time, we find 739 
such halos, including the six main Aquarius Milky Way-like ob- 
jects. We have identi fied substructures in these halos with the 
SUBFIND algorithm dSpringel etal.l200ll) and tracked their or- 
bits by following the position of their most-bound particle. 

We populate these dark matter halos with "galaxies" follow- 
ing a semi-analytic model that uses simple but physically mo- 
tivated laws to track the evol ution of gas cooling, star forma- 
tion a nd feedback processes dLi et alj|2010t [Sta rkenburg et al.l 
120121) . This allows us to derive their baryonic properties such 
as their gas content, stellar mass, etc. Our model simultane- 
ously reproduces the luminosity function, scaling relations and 
chemical content of bri ght as well as dwarf galax ies (for a more 
detailed description see Starkenburg et al. 20l2|). 

When a disk galaxy accretes a low mass companion, it is 
(vertically) heated and puffed up. The increase in the scale 
height AH for a disk of (total, i.e. stellar and gas) mass Md 
and scale length Rd caused by an interaction with a satellite of 
mass M sat may be estimated using analytic arguments to be 
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FIG. 1 . — Left: Gas fraction / gas (blue circles) and galaxy formation efficiency r; ga i (red asterisks) as a function of host mass as predicted by our SA model. The 
dashed blue and solid red curves indicate the respective median trends. Right: The "Spectrum of perturbers" JV(> x) gives the number of encounters with objects 
of a given mass ratio x = M sa t/M v i r (see text for more details). The shaded regions correspond to the 25 and 75 percentiles and were derived from 100 random 
subsamples with each 5 host halos belonging to a given mass range (and are shown for clarity only for the two mass ranges). The black arrows indicate Md/M v i T 
for three representative values of T? ga .i = 5, 30 and 60%, while the thin vertical lines are the SUBFIND mass resolution. 



dToth & Ostrikerll992tlMo et alJ2010h . Here / gas = M gas /M d 
is the gas fraction of the host disk and its inclusion in Eq. (Q~|) 
accounts for the energy that is radiated away and not trans- 
ferred into random motions of disk stars (e.g. iHopkins et al.l 
2008). We have carried out a series of merger experiments on 
the scale of dwarfs (and used analogous simulations of large 
disks bv PVelizauez & White! 1 19991: IVillalobos & Helmil [2008 ; 



iPurcell et al]l200fl iMoster et al.l l2010). and confirm the above 
dependence on the ratio M sat /Md- We have found the propor- 
tionality constant to be a ~ 0.03 when the above expression is 
evaluated at R = 2.5Rd- Below we present two examples of 
such merger simulations and report our results in more detail in 
T. K. Starkenburg et al., in prep. 

Eq. (Q~|) can be re-written in terms of the "disk galaxy effi- 
ciency" of a given halo: 77 gal = M d /(M viT x / bar ), i.e. the 
fraction of baryons collected in the central galaxy compared 
to the total available budget. Here M vu is the virial mass of 
the host halo and /b ar ~ 0.17 is the universal baryon fraction. 
Therefore 

AH a (1 - / gas ) M sat 

/bar 
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Thus three quantities affect the efficiency of disk heating: the 
gas fraction / gas , the galaxy efficiency ry ga i and the mass of the 
perturber compared to that of the host M sat /M v - K . We now 
investigate each of these factors using our models. 

The blue solid circles in the left panel of Fig. [Qshow / gas 
as a function of host halo mass in the SA model. Note that 
the gas content of a galaxy depends strongly on the mass of 
its halo: for objects less massive than 10 10 /i _1 A/© more than 
90% of the baryonic mass assembled onto the central galaxy 
remains as cold gas, revealing how inefficient star formation is 
in (isolated) dwarf galaxies. On the other hand, Milky Way- 
sized objects have typically ~ 10 — 20% of their baryons in 
gas; all these numbers being in reasonably good agreement with 



observations dMcGaugh et al.ll2010b . 

The red asterisks in the left panel of Fig. Q] show ?y ga i as 
function of halo mass. As indicated by the median trend (red 
solid line), halos become increasingly inefficient in collect- 
ing baryons onto galaxies as they become less massive: for 
M vir < lO 10 /*- 1 ^, ?7 ga i ~ 1 - 10%. This is the result 
of a combination of the ef fect of a UV ionizing background 
and of supernova fee dback dLi et al.ll2010t iMaccio et alJ uOlO; 
Okamoto et al. inns. These processes need to be taken into 
account to match the sa tellite luminosity function dGuo et al.l 
|2010t[Moster et al .11201(3) . and explain why dwarf galaxies are 
the most dark matter dominated objects known in the Universe. 

The right panel of Fig.Q]shows the cumulative subhalo mass 
function for our sample of main halos in the Aquarius simula- 
tions for four different ranges of host mass. Since disk heating 
is expected to be more efficient for perturbers that adventure 
close to the center of the host halo, we measure the subhalo 
mass at the first pericenter that is within a distance smaller than 
30% of the virial radius of the host and normalize it to the virial 
mass of the host at that time. The thin vertical lines indicate 
the subhalo resolution, defined by the 20-particle threshold im- 
posed by the SUBFIND algorithm^]. Within the range that is 
well resolved (to the right of the vertical lines), we find that the 
mass spectra of satellites at pericenter are all comparable and 
independent of the virial mass of the host. 

Because the efficiency of galaxy formation 7y ga i depends 
strongly on M v - U (see left panel of Fig. [TJ, at fixed gas content 
the heating produced by satellites is expected to be significantly 
larger for small mass hosts (halos with Af v j r < 10 10 h~ 1 MQ) 
than for Milky Way-like galaxies. To first order this is hinted 

6 Note that for the least massive host halos we are able to resolve fewer 
substructures than for Milky Way-like hosts, and that the SUBFIND algorithm 
is known to underestimate the mass at pericenter, hence the above values are 
lower limits. 
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by the vertical arrows in the right panel of Fig. Q] These ar- 
rows show the mass ratio Md/M v - U for three different values 
of galaxy efficiency: 77 ga i = 5, 30 and 60%, and can be used 
as a guide to determine the number of encounters with satel- 
lites with M sa t ~ Md for a system with a given efficiency. It 
then becomes clear that such encounters are much more com- 
mon for dwarf galaxies, which have lower rj ga i. For example, a 
dwarf galaxy {T] ga \ ~ 5%), experienced on average 1.5 encoun- 
ters with an object of comparable mass in the last 1 1.7 Gyr. On 
the other hand, for a Milky Way-like galaxy whose disk mass is 
~ 10% of the virial value (~ 5% for ri ga \ = 30%), the number 
of encounters with a significant perturber are a factor ~ 15 less 
common. Note that these estimates are somewh at lower than 
derive d in previous work for ~ 10 12 M hosts dPurcell et al.l 
2009) and this could be due to the environment of the Aquarius 
halos. 

3. RESULTS 

Fig. |2]shows, for our model galaxies, AH/ Rd as function of 
host halo mass, normalized to the values expected for a galaxy 
like the Milky Way (with ?7g al = 0.45 and / gas = 0.1) and for 
fixed M sat /M V i r . The red curve indicates the expected change 
when the gas fraction is that of the Milky Way. This shows that 
for a dwarf galaxy populating a 10 9 Ii~ 1 Mq halo, the heating of 
a disk is expected to be ~ 100 times larger than for a galaxy like 
the Milky Way embedded in a 10 12 /i _1 Mq halo. For example, 
even an encounter with a low mass perturber (M sat /M V i r = 
0.05) would be devastating and turn a disky dwarf galaxy into a 
dwarf spheroidal since AH/R d ~ 2.7 for M vir = 1O 9 /i _1 M 
and / gas = 0.1 according to Eq. ©. On the other hand, the 
effect of such an encounter would be nearly negligible in the 
case of a Milky-Way like galaxy. 
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FIG. 2. — Relative increase in the disk's thickness AH/R d as a function 
of the virial mass of the host computed using Eq. {3 for all the main galaxies 
in the Aquarius simulations (black dots). The red curve shows the median ex- 
pected change for / gas = 0.1. All values have been normalized to the heating 
expected for an L„ galaxy, with / gas = 0.1 and ?j ga i = 0.45. The scale on 
the right vertical axis indicates the absolute change in AH/Rd assuming an 
encounter with a M sa t = 0.05M v i r perturber. 



On the scale of Milky Way galaxies the heating is domi- 
nated by subhalos hosting stars, while around smaller hosts 
(M v i r < lO lo ft. _1 M0) the subhalos will generally be dark 
as they fall below the mass threshold imposed by reionization 
and efficient atomic hydrogen cooling to form stars. To con- 
firm that such dark satellites leave imprints on the morpholo- 
gies of dwarf galaxies, we have performed a set of numerical 
experiments. We focus here on two simulations where we var- 
ied the mass ratio between the disk and the satellite, but took 
M sat /M v i r = 0.2 comparable to what has been used in pre- 
vious wo rk ( Kazantzidis et al]|200aiVillalobos & Helmill2008l: 
iPurcell et aU2009tlMoster et alj|2010h . T he satellite follows an 
NFW profile with concentration c = 18.7 dMufloz-Cuartas et alj 
l201lb . Our disk galaxies are purely stellar, they have Md = 

0. 008 and 0.04 x M V \ T , and are embedded in a Hernquist halo 
with mass M vn = 10 10 /i -1 M Q , and scale-radius a = 9.3 kpc, 

1. e. ?7g a i ~ 5% and 23% respectively. The disks are radially 
exponential with scale-length Rd = 0.67ft- -1 kpc, and verti- 
cally they follow sech 2 (z/2,2ro), wit h zp = 0.05Rd. T he inter- 
nal kinematics are set-up following Hernquist (1993), and the 
disks are stable (with Toomre parameters Q > 1). 

We put the satellite on a fairly radial orbit with r apo /r peT i — 
40, starting from a distance of ~ 23 kpc, and found that it is 
completely disrupted after three close passages, i.e. in ~ 1.5 
Gyr. Fig. [3] shows the final surface brightness profiles of the 
heavy and light disks in the top and bottom panels respectively. 
This figure evidences that significant heating has taken place 
and even led to important changes in the morphology of the host 
galaxy. This is expected since although we simulated minor 
mergers in terms of virial mass ratios, these are major mergers 
from the perspective of the dwarf galaxy, as M sat /Md = 5 and 
25 respectively. 

In the case of gas-rich systems -which are a majority at the 
low mass end, encounters with dark satellites will be less effi- 
cient at changing the structure of the host dwarf galaxy, because 
much of the orbital energy will be absorbed by the gas, lead- 
ing to less vertical heating. However, we may expect that such 
encounters may induce star formation events, and thus, albeit 
indirectly, lead to significant changes in the characteristics of 
these galaxies (T. K. Starkenburg et al., in prep). 

To establish whether observations support that disks of dwarf 
galaxies are thicker than those of larger systems, we have com- 
piled measurements of the thickness of stellar disks (quanti- 
fied by the apparent axis ratio, b/a) for a wide range of galaxy 
masses. Although the observed b/a is not a measurement of 
the intrinsic shape of the disk, if one assumes random orienta- 
tions on the sky, the two are directly related. In our literature 
search we have carefully selected isolated late-type galaxies to 
avoid any morphology-luminosity trend that ma y be driven b y 
environmental interactions (such as discussed in Maver 2010). 

The top panel of Figure|4]shows the distribution of optical (r 
or R band) b/a as a function of circular velocity (V c = Wso/2). 
The latter provides a measure of the dynamical mass of the 
galaxy and its dark matter halo. We plot here data fo r two 
galaxy samples: HOPCAT (black dots: lDovle et al.l l2005) con- 
taining the optical counterparts of ~ 3600 HIPASS sources, and 
for a set of isolated nea rby late type galaxies (blue asterisks; 
iKarachentsev et af]|2004 tidal index < 0, RC3 morpholog- 
ical type > 0). The magenta square shows the median valu e 
for a subsample of the 101 dwarf galaxies dGeha et al.l l2006). 
where we have selected only those objects with no companions 
within 1 Mpc projected distance and (g — r) < 0.55. The bot- 
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FIG. 3. — The left panels show the initial surface brightness profiles for two of our simulated dwarf galaxies with M,j/Afvir = 0.008 and 0.04 (bottom and top 
panels, respectively). The panels on the right correspond to the final stellar distributions after these disks merged with a dark satellite of mass M sa t = 0.2M v ; r , 
and are shown after 6 Gyr of evolution (i.e well after the merger has taken place, so the system appears to be relaxed again). 
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FIG. 4. — Apparent optical axis ratios b/a for various samples of isolated, 
late-type galaxies from the literature. Top panel: HOPCAT (Doyle et al. 2005; 
black dots), a sample of nearby galaxies selected from Karachentsev et al. 
(2004; blue asterisks) and the median of the sample of dwarf galaxies pre- 
sented by Geha et al. (2006; magenta square with vertical bar showing the 
25-75% percentiles). Bottom panel: Edge-on galaxies from Yoachim & Dal- 
canton (2006). These are shown separately because the b/a corresponds to 
their true axis ratio, i.e. zq/ R d . The lines indicate, for each sample, the me- 
dian 6/ a at a given circular velocity V c . 



torn panel of this Figure shows the intrinsic thickness defined as 
the rat io of scale-height to scale-le ngth for a sample of edge-on 
disks (I Yoachim & Dalcantonll2006l) . 

The median b/a trends with galaxy circular velocity are in- 
dicated separately f or each sample by the bla ck solid (HOP- 
CAT), blue dashed (iKarachentsev et al.l 12004) and red dotted 
dYoachim & Dalcantonl 20061) curves. Each set clearly shows 
that the axis ratios increase with decreasing circular velocities. 
In other words, stellar disks become thicker as we move to- 
wards less massive galaxies, in qualitative agreement with our 
expectations based on the analysi s of disk heating by substruc - 
ture on isolated galaxies (see also lSanchez-Janssen et al.l20 10). 

4. DISCUSSION 

We have demonstrated that the dynamical effects of dark 
satellites on disky dwarf galaxies are much more dramatic than 
on galaxies like the Milky Way. Mergers with M sat > Md 
are not very common for z < 3 but they occur much more 
frequently than on the L»-galaxies scale. As an example, we 
have simulated a merger with M sat /M v i r = 0.2 for a dwarf 
with Md = 8x 10 7 /i _1 M Q in stars, i.e. slightly more massive 
than the Fornax dwarf galaxy, and found that its morphology 
changed from disky to spheroidal. This might be a plausible 
path for the formation of dSph systems in isolation (if the dwarf 
was gas poor, which is rare in our models but not unlikely). 
This channel might also be relevant for the dSph satellites of 
our Galaxy, provided such encounters would have taken place 
just before the system fell onto the potential well of the Milky 
Way (since further gas accretion would thus be prevented). 

Most of the galaxies on the scales of dwarfs are, however, 
gas-rich. In that case, encounters with dark satellites can trig- 
ger starbursts, which might explain the presence of seemingly 
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isolated dwarfs undergoing major star formation events with- 
out an apparent trigger. Depending on the characteristics of the 
encounter, such starbursts will vary in amplitude. We are cur- 
rently performing hydrodynamical simulations to characterize 
this process (T. K. Starkenburg et al. in prep). 

Additionally, other processes exist that can influence the 
morphologies of dwarf galaxies. For example, binary mergers 
between disky dwarfs can result in the formation of spheroidal 
systems dKazantzidis et al.l201 lb . although such events are rare 
(see Fig. 1). On the other hand, on the scale of dwarf galax- 
ies physical processes affecting gas may also lead to thicker 
systems. For example, the presence of a temperature floor in 
the interstellar medium at T ~ 1Q A K introduced by e.g. a 
UV background, implies that gas pressure becomes compara- 
ble to rotational support for small dark matter halos. Stars 
formed in such systems wo uld thus be born i n puffi er config- 
urati ons as demonstrated by iKaufmann et al.l (|2007) (see also, 
e.g. iRobertson & Kravtsovll2008l) . 

Yet, we have shown here that a distinctive imprint on dwarf 
galaxies will be left by dark satellites in the context of the 
ACDM cosmological paradigm. Such dark satellites are ex- 
pected to make the stellar disks of isolated dwarf galaxies sig- 
nificantly thicker than those of ~ L„ galaxies. We have indeed 
detected such a trend on three different observational samples 
of isolated late-type galaxies on the nearby Universe. We may 
have identified a new mechanism to explain the morphologies 
of dwarf galaxies. 
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